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1.0  Introduction 


This  annual  report  summarizes  research  for  the  period 
15  April  1979  to  14  March  1981  on  the  measurement  and  analy¬ 
sis  of  dense  non-ideal  plasmas  generated  by  a  compact,  ex¬ 
plosively  driven  device.  In  previous  programs  dating  back 
to  1976,  Artec  Associates  has  pursued  a  program  of  basic 
research  in  pulsed  plasma  MHD  (References  1,  2  and  3).  As 
a  result,  our  knowledge  of  the  properties  of  dense  non-ideal 
plasmas  has  been  greatly  expanded.  We  have  generated  plasma 
states  intermediate  to  those  achieved  by  chemical  and  nuclear 
explosive  means;  we  have  experimentally  measured  plasma  elec¬ 
trical  conductivity  and  flow  velocity  histories;  and  we  have 
helped  define  the  theoretical  framework  for  predicting  con¬ 
ductivities  in  the  dense  non-ideal  regime. 

The  results  of  the  present  work  have  extended  our  know¬ 
ledge  and  now  form  the  underlying  theoretical  and  experimen¬ 
tal  basis  for  pulsed  plasma  MHD  as  a  compact  primary  power 
source  capable  of  megajoule  electrical  pulses  at  gigawatt 
power  levels  and  extremely  high  power  densities.  We  have 
satisfactorily  characterized  the  plasma  flow  and  have 
achieved  good  agreement  between  measured  and  calculated 
plasma  conductivities.  The  research  work  summarized  in 
this  report  includes  the  results  of  our  last  three  plasma 
experiments  (Shots  130-5,  130-6  and  130-7),  as  well  as  the 
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results  of  i  comprehensive  analysis  of  all  six  plasma  ex¬ 
periments  that  have  used  the  high  efficiency  plasma  source 
(Shots  130-2  through  130-7).  The  experimentally  measured 
flow  velocities,  conductivities  and  other  diagnostics  are 
analyzed  and  compared  to  those  predicted  from  hydrodynamic 
calculations  of  the  plasma  generation  and  subsequent  channel 
flow. 


The  plasma  source  used  in  these  latest  experiments 
drives  a  30  to  40  diameter  column  of  high  velocity,  high 
conductivity  plasma  with  an  average  hydrodynamic  power  on 
the  order  of  a  100  gigawatts.  Passing  the  plasma  through 
a  magnetic  field  of  strong  permanent  magnets  we  have  ex¬ 
perimentally  achieved  peak  MHD  power  densities  of  3.6  tera- 
watts/m3-tesla2 . * 

The  objective  of  our  next  series  of  experiments  is  to 
demonstrate  peak  electrical  power  in  the  gigawatt  range 
using  an  externally  generated  10  tesla  field.  The  success¬ 
ful  completion  of  these  experiments  will  demonstrate  an 
electrical  output  of  3.5  kilojoules  per  kilogram  of  system 
weight  in  an  expendable  single  pulse  system.  Further  im¬ 
provements  in  the  output  of  the  plasma  source  and  configura¬ 
tion  of  the  explosives  can  double  or  triple  this  energy  out¬ 
put  per  system  weight.  Further  increases  in  magnetic  field 


*  power  density  in  terms  of  the  effective  volume  between  the 
electrodes  and  the  average  magnetic  field  strength  in  this 
volume. 
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strength,  for  example  by  self-excitation,  can  yield  even 
larger  electrical  power  and  energy  outputs  until  extrac¬ 
tion  rates  begin  to  seriously  interact  with  and  brake  the 
plasma  flow. 

We  have  considered  a  repetitively  pulsed  cartridge 
loaded  configuration  and  estimate  an  electrical  output  of 
50  joules  per  kilogram  of  system  weight  per  pulse  for  a 
completely  contained  system.  The  repetition  rate  of  such 
a  contained  system  would  be  commensurate  with  rapid  fire 
gun  technology  and  would  be  dependent  on  cartridge  size. 

The  high  power  density  characteristics  of  pulsed 
plasma  MHD  lend  themselves  to  a  number  of  military,  re¬ 
search  and  commercial  applications.  Military  applications 
include  use  as  a  prime  power  source  in  the  externally  or 
partially  self-excited  mode  for  directed  energy  weapons 
and  electromagnetic  guns.  By  virtue  of  their  compact  size 
and  high  energy  and  power  per  unit  system  weight,  these 
generators  are  being  considered  as  a  lightweight  power 
supply  for  space  based  weapons,  surveillance  and  various 
forms  of  electronic  warfare.  These  generators  may  also 
be  suitable  for  strategic  communications  from  low  earth 
orbit  to  submerged  submarines. 

The  pulsed  plasma  MHD  generator  can  provide  very  high 
magnetic  fields  (in  the  self-excited  mode)  or  act  as  a 
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current  source  for  flux  compressors  in  various  research 
and  testing  applications  including  electromagnetic  guns 
for  hypervelocity  research.  On  site  pulse  testing  of 
large  power  transformers  and  field  resistance  welding 
of  large  structures  have  been  identified  as  promising 
commercial  applications. 


2.0  Plasma  State 


In  our  experiments  the  working  gas  is  energized  at 
high  pressures  (3  GPa  to  100  GPa)  in  a  compact  explosively 
driven  device  that  efficiently  converts  the  chemical  energy 
of  the  explosive  to  the  working  gas/plasma.  Conversion 
efficiencies  in  excess  of  20%  have  been  achieved.  The 
plasma  is  then  expanded  down  a  long  channel  to  form  a  high 
velocity  extended  duration  pulse.  The  length  of  the  plasma 
column  is  typically  30  to  40  channel  diameters  and  the  peak 
flow  velocity  is  typically  25  to  30  km/s. 

The  argon  plasma  states  in  which  we  have  conducted 
most  of  our  experiments  are  shown  in  Figure  1  as  contours 
of  electrical  conductivity  on  an  internal  energy  versus 
pressure  diagram.  The  conductivities  are  calculated  using 
the  Rogov  formulation  (Reference  4)  based  on  a  modified 
Debye-Huckel  equation-of-state  description.  These  calcu¬ 
lations  are  described  in  further  detail  in  References  2 
and  5.  The  hydro  code  flow  calculations  of  our  plasma 
experiments  used  the  1980  Livermore  argon  equation-of- 
state  (Reference  6) .  In  the  range  of  conditions  encoun¬ 
tered  in  the  channel  flow,  the  Debye-Huckel  equation-of- 
state  yields  state  parameters  within  a  few  percent  of  the 
Livermore  equation-of-state.  Thus  the  conductivities  de¬ 
rived  from  the  flow  calculations  using  the  Livermore 
equation-of-state  are  consistent. 
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Figure  1.  Argon  Plasma  States  and  Electrical 
Conductivity  (kS/m) 


The  peak  plasma  conditions  generated  by  the  high 
efficiency  plasma  source  are  characterized  by  100  GPa 
pressures  and  near  metallic  densities.  In  this  regime 
the  Debye-Huckel  equation-of-state  is  not  applicable 
and  we  have  used  only  the  Livermore  equation-of-state 
to  calculate  plasma  source  dynamics. 

In  all  of  our  experiments  using  various  working 
gases  and  plasma  source  geometries,  the  bulk  electrical 
conductivities  measured  or  calculated  in  the  plasma 
channel  range  from  50  kS/m  at  peak  flow  velocities  of 
30  km/s  to  15  kS/m  at  the  end  of  a  highly  expanded  pulse. 
The  corresponding  pressure  range  is  1  GPa  to  50  MPa.  A 
solid  theoretical  foundation  is  essential  for  generating 
and  measuring  plasmas  in  this  dense,  non-ideal  regime. 

We  are  confident  of  the  1980  Livermore  equation-of-state 
for  argon  and  the  Rogov  formulation  of  electrical  con¬ 
ductivity  and  we  are  consistently  getting  good  agreement 
between  measured  and  calculated  flow  parameters. 

Table  1  summarizes  typical  plasma  states  encountered 
in  the  plasma  source  and  plasma  channel.  These  estimates 
are  based  on  the  calculations  and  experiments  described 
in  Section  4.0. 


15 


Pressure  Energy  Density  Conductivity 

(GPa) _ (MJ/kg)  (Mq/m3) _ (kS/m) 


Plasma  Source  and  Plasma  Channel 


3.0  Recent  Plasma  Experiments 


3.1  Description  of  the  Experiments 

The  high  efficiency  plasma  source  is  shown  in  cross- 
section  in  Figure  2.  Octol  explosive  is  cast  around  a 
steel  lined  annular  region  containing  a  pressurized  working 
gas.  By  means  of  the  wave  shaper  the  detonation  front  is 
caused  to  pass  simultaneously  along  the  outer  and  inner 
liners  of  the  annular  region.  The  impact  of  the  liners 
forms  a  dynamic  seal  which  progresses  along  the  annulus 
at  the  explosive  detonation  velocity  driving  a  strong, 
high  pressure  shock  into  the  working  gas.  The  result  is 
an  efficient  conversion  of  explosive  energy  into  the 
strongly  shocked  working  gas/plasma. 

The  plasma  converges  into  the  hemispherical  end 
cavity  and  as  it  stagnates  it  is  further  energized  by 
the  central  piston  driven  by  the  inner  core  of  explosive. 

It  thereupon  expands  into  the  partially  evacuated  plasma 
channel . 

Air  was  used  as  the  working  gas  on  Shot  130-5  to 
provide  a  complement  to  the  previous  work  using  xenon 
and  argon.  Argon  was  used  as  the  working  gas  on  Shots 
130-6  and  130-7,  and  the  end  cavities  were  modified  as 
shown  in  Figure  3.  Calculations  have  shown  that  the 
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Figure  2.  Cross-Section  of  Explosive  Plasma  Source 


Shot  ‘ 30-6 


Figure  3.  Modified  End  Cavity  Sections  of 
Explosive  Plasma  Sources 


100  gigapascal  pressures  in  the  convergence  section  lead 
to  early  closure  of  the  channel  inlet.  To  ensure  that 
most  of  the  plasma  is  driven  into  the  channel,  the  last 
two  experiments  (Shots  130-6  and  7)  were  designed  with 
a  short  nozzle  inlet  connecting  the  cavity  to  the  channel. 
Active  tamping  of  the  cavity  was  provided  by  additional 
explosive  as  shown  in  Figure  3 . 

Three  types  of  diagnostics  were  used  to  measure  the 
plasma  properties  in  the  channel  as  illustrated  in  Figure 
4.  The  plasma  flow  velocity  history  was  monitored  by 
means  of  an  open  circuit  Faraday  generator.  When  a  plasma 
passes  through  a  transverse  magnetic  field,  a  voltage  is 
produced,  when  the  flow  velocity,  magnetic  field  and  pick¬ 
up  electrodes  are  mutually  orthogonal,  the  open  circuit 
voltage  is  given  by 

vQ  =  Bub  (la) 

where  b  is  the  interelectrode  distance  in  meters,  u  is  the 
velocity  in  m/s,  and  B  is  the  magnetic  field  strength  in 
teslas.  With  a  known  magnetic  field  and  known  electrode 
spacing  the  measured  voltage  provides  the  plasma  velocity 
history  by 


u  (t ) 


Vo  (t) 
Bb 


(lb) 
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Figure  4.  Typical  Plasma  Channel 
Diagnostic  Station 


In  a  similar  manner  the  plasma  resistance  history, 
and  hence  effective  plasma  conductivity,  was  monitored 
by  means  of  a  loaded  Faraday  generator.  In  this  case, 
the  generator  is  provided  with  a  10  mft  resistive  load 
which  roughly  matches  the  plasma  resistance.  The  vol¬ 
tage  appearing  across  the  load  electrodes  is  given  by 


vT  =  v  -  R  I  -  L  4r 
L  o  P  p  dt 


(2) 


where  I  is  the  load  circuit  current  and  Rp  is  the  plasma 
resistance.  The  plasma  resistance  history  is  obtained 
by  measuring  vL  vQ ,  and  I.  VQ  is  obtained  from  the  ve¬ 
locity  gage,  I  is  obtained  from  a  search  coil  inductively 
coupled  to  the  load  and  Lp  is  estimated  from  ^  at  t=0. 


The  effective  plasma  conductivity  is  determined  from 
the  plasma  resistance  by 


where  r  is  an  electrode  geometry  factor.  This  was  measured 
in  the  laboratory  with  an  electrode  and  channel  mock-up  by 
passing  a  suitable  current  through  an  electrolyte  of  known 
conductivity. 


The  plasma  conductivity  history  was  also  measured 
with  pickup  coils  placed  near  the  magnet  ends  as  shown. 
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When  the  plasma  flows  through  an  axial  gradient  in  the 
B-field,  there  is  an  axial  gradient  in  the  induced  Fara¬ 
day  voltage,  and  this  results  in  circulating  currents  in 
the  plasma.  A  search  coil  placed  outside  the  channel 
will  have  a  voltage  induced  in  it  by  the  time  variations 
in  the  magnetic  fields  produced  by  these  eddy  currents. 

The  outputs  of  the  search  coils  are  electronically  inte¬ 
grated,  so  they  directly  measure  the  field  perturbations 
caused  by  the  eddy  currents: 

v  =  NAAB/t  (4) 

where  A  is  coil  area,  N  is  number  of  turns,  and  t  is  the 
integrator  time  constant.  The  magnetic  field  perturbation 
AB  is  a  function  of  magnetic  Reynolds  number,  being  theore¬ 
tically  linear  with  Rm  at  low  values.  Thus 

AB  =  f(Rm)  =  f(p0aub)  (5) 

The  search  coils  are  calibrated  in  the  laboratory  by 
propelling  a  metal  bar  of  known  conductivity  through  a  pre¬ 
cise  mockup  of  the  shot  geometry  to  obtain  coil  output  as 
a  function  of  magnetic  Reynolds  number.  This,  together  with 
the  measured  velocity  history,  is  used  to  obtain  the  conduc¬ 
tivity  history  from  the  measured  coil  output  on  each  shot. 
Further  details  on  the  design  and  calibration  of  the  plasma 
diagnostics  are  given  in  Reference  2. 


3.2  Experimental  Results 


The  first  three  plasma  experiments  which  used  the 
high  efficiency  plasma  source  have  been  reported  in 
Reference  3.  These  are  Shot  130-2  with  xenon  as  the 
working  gas;  Shot  130-3  with  argon;  and  Shot  130-4  which 
used  a  scaled-up  version  of  the  130-3  plasma  source. 

Shot  130-1  which  used  a  different,  less  energetic  plasma 
source  is  not  included  in  the  present  analysis.  The 
following  sections  present  for  the  first  time  the  re¬ 
sults  of  Shots  130-5,  6  and  7. 

The  measured  flow  velocity  history,  conductivity 
histories  and  normalized  load  power  history  for  130-5 
(with  air  as  the  working  gas)  are  shown  in  Figures  5, 

6  and  7.  The  velocity  history  suaqests  a  somewhat  less 
energetic  plasma  than  observed  in  previous  experiments. 

The  conductivity  is  initially  comparable  to  that  of  pre¬ 
vious  experiments  but  decays  more  rapidly  with  time  and 
with  distance  down  the  channel. 

The  measured  load  power  history,  like  the  flow  ve¬ 
locity  history,  is  lower  than  previous  experiments.  Based 
on  an  analysis  of  the  plasma  source  data,  we  believe  that 
the  mass  of  plasma  driven  into  the  channel  was  less  than  in 
previous  shots.  The  analysis  in  part  based  on  the  cavity 
X-ray  (Figure  8)  concludes  that  because  of  the  higher  shock 
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Figure  8.  Digitized  Representation  of  the 
120  MEV  Flash  X-Ray  of  the  130-5 
Plasma  Source  Cavity 


compressibility  of  air  the  plasma  source  cavity  volume 
collapses  prematurely  and  traps  a  substantial  mass  of 
plasma  isolating  it  from  the  channel  inlet.  The  plasma 
densities  in  the  channel  flow  are  therefore  reduced  and 
consequently  the  available  hydrodynamic  power  is  reduced. 

The  results  for  Shot  130-6  (with  argon)  are  given 
in  Figures  9,  10  and  11.  Analysis  of  the  cavity  X-ray 
(Figure  12)  again  reveals  that  a  substantial  fraction  of 
the  plasma  is  trapped  by  a  premature  collision  of  the 
central  piston  with  the  cavity  walls.  The  additional 
volume  of  the  nozzle  inlet  results  in  lower  cavity  pres¬ 
sures  and  the  central  piston  motion  is  not  reversed  in 
time.  The  plasma  density  in  the  subsequent  channel  flow 
is  thus  reduced  and  as  a  result  conductivity  and  flow 
power  levels  are  low. 

The  cavity  for  Shot  130-7  (with  argon)  was  rede¬ 
signed  to  force  most  of  the  plasma  into  the  channel  by 
eliminating  the  regions  in  the  cavity  where  plasma  can 
be  trapped.  The  central  piston  is  allowed  to  continue 
into  the  tapered  channel  inlet.  This  design  is  expected 
to  result  in  a  somewhat  lower  energy-entropy  state  for 
the  plasma.  The  cavity  X-ray  (Fiqure  13)  shows  that  the 
cavity  collapsed  in  the  manner  expected. 
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Shot  130-6  (Argon) 
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3.4  ys  After  First  Motion  of  the 
Front  of  the  Central  Piston 


Figure  12. 


Digitized  Representation  of  the 
120  MEV  Flash  X-Ray  of  the  130-6 
Plasma  Source  Cavity 


6.9  us  After  First  Motion  of  the 
Front  of  the  Central  Piston 


Figure  13 


Digitized  Representation  of  the 
120  MEV  Flash  X-Ray  of  the  130-7 
Plasma  Source  Cavity 


The  results  of  the  experiment  are  shown  in  Figures 
14,  15  and  16.  The  general  leve3  of  the  velocity  history 


is  down  from  that  of  130-2  or  130-3  in  accordance  with 
the  expectation  of  a  lower  energy-entropy  reservoir  state. 
However,  conductivity  levels  are  high  as  would  be  expected 
with  a  good  fraction  of  the  source  plasma  driven  into  the 
channel.  The  effective  conductivity  is  measured  by  the 
load  electrodes  and  as  in  all  previous  experiments  is 
lower  than  would  be  measured  by  the  search  coil  conducti¬ 
vity  gages.  The  normalized  load  power  levels  are  somewhat 
higher  than  those  of  Shot  130-3,  independently  confirming 
the  high  levels  of  measured  conductivity. 
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Figure  16.  Normalized  Load  Power  History 
For  Shot  130-7  (Argon) 
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4.0  Analysis  of  Plasma  Experiments 

Upon  completion  of  the  three  plasma  experiments  de¬ 
scribed  in  the  previous  section,  we  made  a  comprehensive 
analysis  and  comparison  of  the  plasma  measurements  from 
all  six  shots  (130-2  through  130-7)  that  used  the  high 
efficiency  plasma  source.  To  complement  this  analysis 
we  developed  a  computational  model  for  the  plasma  source 
and  channel  hydrodynamics  and  used  the  Rogov  formulation 
for  conductivity  to  estimate  conductivity  histories.  A 
summary  of  the  major  parameters  of  the  six  plasma  experi¬ 
ments  is  given  in  Table  2. 

4.1  Experimental  Comparisons 

Three  of  the  experiments  (130-2,  3  and  5)  compare 
three  different  working  gases- -xenon ,  argon  and  air-- 
using  otherwise  identical  plasma  sources,  plasma  channels 
and  diagnostics.  The  loading  density  of  each  experiment 
was  adjusted  to  provide  the  same  mass  of  plasma  at  approx¬ 
imately  the  same  pressures.  During  the  operation  of  the 
plasma  source,  the  xenon  and  argon  exhibit  about  the  same 
compressibility.  The  compressibility  of  the  air  at  condi¬ 
tions  prevailing  in  the  plasma  source  is  significantly 
higher  and  from  analysis  of  the  plasma  source  diagnostics 
we  believe  that  the  mass  of  plasma  driven  down  the  channel 
was  considerably  less  than  in  the  xenon  or  argon  experiments 
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This  conclusion  is  further  strengthened  by  the  lower  flow 
velocity  and  normalized  load  power  observed  in  the  plasma 
channel.  Otherwise  these  three  experiments  represent  a 
comparison  of  the  conductivities  of  the  three  working  gases. 

The  measured  flow  velocity  histories,  conductivity 
histories  and  normalized  load  power  histories  for  the 
three  shots  are  shown  in  Figures  17,  18,  19  and  20.  Xenon 
conductivities  average  about  30%  higher  than  those  of 
argon.  Based  on  our  best  estimates  and  taking  account  of 
plasma  source  performance,  it  appears  that  air  conductivity 
may  be  as  high  or  higher  than  that  of  argon  at  the  same 
thermodynamic  conditions. 

The  channel  shock  velocities  for  all  six  experiments 
are  listed  in  Table  2.  From  previous  work  (Reference  7) 
we  have  noted  that  framing  camera  records  of  the  flow  have 
often  shown  a  highly  canted  non-uniform  shock  front.  At 
least  the  first  few  diameters  of  channel  flow  is  not  one¬ 
dimensional  and  for  this  reason  we  have  made  no  attempt 
to  relate  the  peak  measured  flow  velocity  to  the  one-di¬ 
mensional  Hugoniot  flow  velocity  for  the  observed  shock 
velocity . 

In  130-4  the  dimensions  of  the  plasma  source  were  scaled 
by  a  factor  of  1.8  over  those  of  130-3.  The  channel  and 
diagnostics  were  unchanged.  The  objective  of  130-4  was 
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Figure  19.  Conductivity  Histories  For 
Shots  130-2,  3  and  5 
(Conductivity  Gage  2) 


to  extend  the  duration  of  channel  flow  by  increasing  the 
duration  of  the  high  channel  inlet  conditions.  The  ve¬ 
locity,  conductivity  and  normalized  load  power  histories 
for  the  two  experiments  are  shown  in  Figures  21,  22  and 
23  and  it  is  evident  that  performance  gains  with  the 
larger  source  are  modest.  There  are  two  leading  hypo¬ 
theses  to  account  for  this. 

1.  The  channel  inlet  is  collapsed  by  the  85 
GPa  plasma  source  pressures  thereby  cutting  off  the 
plasma  flow  at  approximately  the  same  time  in  both  ex¬ 
periments.  This  effect  is  predicted  and  roughly  quanti¬ 
fied  by  the  2-D  code  calculations  supporting  the  Lawrence 
Livermore  Voitenko  compressor  experiments  (Reference  8) 
which  gave  peak  channel  inlet  conditions  very  close  to 
our  experiments. 

2 .  Boundary  layer  growth  in  the  high  pressure 
region  of  the  channel  near  the  plasma  source  is  driven 
by  wall  ablation  and  closes  within  10  to  20  diameters 

of  the  shock  front.  We  are  currently  incorporating  boun¬ 
dary  layer  effects  into  our  flow  calculations  to  deter¬ 
mine  whether  the  flow  conditions  just  downstream  of  the 
plasma  source  can  generate  a  sufficiently  blown  boundary 
layer  to  affect  the  flow  in  the  ways  we  observe. 
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21.  Flow  Velocity  Histories  For 
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Figure  22.  Conductivity  Histories 
For  Shots  130-3  and  4 


The  plasma  source  used  in  the  130-7  experiment  rep¬ 
resents  an  improvement  in  the  design  of  the  end  cavity 
and  channel  inlet.  The  comparison  of  channel  diagnostics 
for  130-3  and  130-7  (Figures  24,  25  and  26)  illustrates 
the  effect  of  these  changes.  As  judged  by  conductivity 
and  load  power  histories  the  improvements  in  channel  flow 
conditions  are  modest  even  though  an  analysis  of  the  plas¬ 
ma  source  cavity  collapse  would  indicate  that  considerably 
more  plasma  has  to  be  driven  into  the  channel.  The  rela¬ 
tive  insensitivity  of  pulse  duration  to  improvements 
in  plasma  source  design  may  be  another  indication  of  channel 
flow  controlled  by  blown  boundary  layer  effects. 

As  reported  in  our  previous  work  (References  1,  2  and 
3)  the  conductivity  measured  by  the  search  coil  conductivity 
gages  is  invariably  higher  than  the  effective  plasma  con¬ 
ductivity  calculated  from  the  load  electrode  current  and 
voltage  traces.  The  conductivity  gage  measures  the  magnetic 
flux  from  eddy  currents  within  the  plasma  and  therefore 
tends  to  weight  the  regions  of  higher  conductivity.  The 
effective  conductivity  derived  from  load  electrode  measure¬ 
ments,  however,  would  tend  to  weight  the  more  resistive 
boundary  layer  regions  and  be  responsive  to  the  details  of 
electrical  contact  between  plasma  and  electrodes.  We  con¬ 
tinued  to  observe  this  difference  in  Shots  130-5  through 
130-7.  We  use  the  lower  and'more  relevant  effective  con¬ 
ductivity  values  in  our  estimates  of  MHD  power  generation. 
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Figure  25.  Conductivity  Histories 
For  Shots  130-3  and  7 


Figure  26.  Normalized  Load  Power  Histories 
For  Shots  130-3  and  7 
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The  normalized  load  power  curves  such  as  shown  in 
Figures  23  and  26  are  derived  by  multiplying  the  measured 
load  current  by  the  measured  voltage  then  dividing  by  the 
square  of  the  measured  magnetic  field  strength.  The  gen¬ 
eration  volume  as  defined  by  the  plasma  volume  accessible 
to  the  load  electrodes  is  2.3  cm3.  For  Shot  130-7  the 
peak  normalized  power  density  is  3.6  TW/m*-T-  (Figure  16) 
and  the  normalized  power  density  averaged  over  the  first 
2  5  us  of  flow  is  2.4  TW/nr-T2  . 

While  this  power  density  was  measured  for  a  field 
strength  of  about  0.1  tesla  we  expect  that  it  will  be  rea¬ 
sonably  constant  up  to  field  strengths  of  at  least  10  tesla. 
Much  of  this  confidence  is  based  on  the  high  currents  and 
field  amplifications  observed  in  previous  work  on  self-exci- 
-  -.fir.  I  ?o  f  rencos  2,  2  and  9).  Current  in  the  self-excited 
■  :  increased  as  expected  until  the  process  quenched  as  a 

r-'  suit  :  the  limited  pulse  length. 

Consider  a  resistive  and  inductive  load  circuit  with 
the  current  given  by  the  relation 

(Ll  +  Lp)  C^+Rp)I  -  Bub  =  0  (6) 

and  the  lend  •.’•'Itam:  bv 


=  L  ^ 

L  dt 


V 


(7) 
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For  load  and  plasma  inductances  typical  of  our  experi¬ 
ments,  the  inductive  voltage  drop  can  be  neglected  after 
the  first  few  diameters  of  flow.  Thus,  for  a  matched 
load,  the  load  power  is 


PL  -  V 


(Bub) 2 
4RL 


(8a) 


Expressing  the  plasma  resistance  in  terms  of  conductivity, 
the  load  power  becomes 


B2u2bAa 
*L  4 

B  2  uARm 

where  Rra  =  u0aub  (magnetic  Reynolds  number) 
and  A  =  effective  electrode  area 


(8b) 


For  a  10  tesla  field,  large  electrodes  and  an  average 
Rjn  of  10  (typical  of  our  experiments)  ,  the  average  elec¬ 
trical  power  extracted  would  be  10  gigawatts  with  a  peak 
of  nearly  twice  this  value.  The  average  electrical  power 
extracted  is  about  10%  of  the  average  hydrodynamic  power 
available  in  the  flow  and  might  be  expected  to  cause 
moderate  flow  perturbations. 


55 


4.2  Calculation  of  Channel  Flow  and  Conductivity 

Histories 

4.2.1  Calculational  Method 

We  have  utilized  an  advanced  1^-D  hydro  code 
with  interactive  wall  motion  to  model  the  energizing  of 
the  plasma  and  subsequent  channel  flow.  The  collapse 
of  the  plasma  source  cavity  is  made  to  conform  closely 
to  a  number  of  experimentally  observed  constraints.  This 
approach  is  valid  because  of  the  high  sound  speeds  charac¬ 
teristic  of  the  plasma  in  the  source.  It  gives  a  good 
account  of  overall  flow  dynamics  and  allows  a  rapid  eval¬ 
uation  of  parameters  controlling  the  flow. 

Our  calculations  are  based  on  the  1-D  streamtube 
finite  difference  equations  of  Reference  10.  The  working 
gas  is  modeled  by  a  28  term  polynomial  fit  to  the  1980 
Livermore  argon  equation-of-state  (Reference  6)  and  is  the 
best  available  data  for  the  entire  range  of  plasma  states 
encountered  in  the  plasma  source  and  subsequent  channel 
flow.  The  channel  is  modeled  by  a  series  of  wallpoints 
(Figure  27)  whose  radial  motion  is  controlled  by  the  local 
channel  pressure  and  by  the  shock  Hugoniot  characteristics 
of  the  material  associated  with  the  wallpoints.  Channel 
wallpoints  at  diagnostic  locations  corresponding  to  those 
of  the  experiments  are  used  to  collect  the  appropriate 
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hydrodynamic  variable  data  for  direct  comparison.  In  the 
case  of  conductivity  gage  locations,  the  pressure  and  in¬ 
ternal  energy  histories  are  saved  and  used  in  a  post-pro¬ 
cessor  code  to  compute  a  conductivity  history  based  on  a 
table  lookup  compiled  from  the  Rogov  formulation  of  plasma 
conductivity  (References  4  and  5).  Load  voltage  and  cur¬ 
rent  are  also  computed  by  the  post-processor  code  from  a 
saved  velocity  history  and  derived  conductivity  history. 
These  histories  are  used  in  a  numerical  solution  of  the 
load  circuit  equation  (6)  to  derive  current  and  voltage 
histories . 


The  air  in  the  channel  ahead  of  the  working  gas 
is  modeled  by  an  air  Hugoniot  boundary  condition. 

Modeling  the  operation  of  the  plasma  source  is 
the  key  to  achieving  a  reasonable  channel  flow  calculation. 
Referring  to  Figure  27,  the  calculation  is  begun  when  the 
central  piston  first  begins  to  move.  The  state  of  the 
working  gas  in  the  cavity  is  determined  fully  by  specifying 
two  state  variables.  The  density  of  the  gas  is  calculated 
from  the  volume  of  the  cavity  and  the  estimated  mass  of 
gas  in  the  cavity.  The  entropy  of  the  gas  is  assumed  to 
be  that  of  the  driver  gas  when  brought  to  rest  by  a  strong 
shock  reflection.  The  initial  cavity  state  is  then  taken 
to  lie  on  the  isentrope  passing  through  this  reflected 
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shock  state.  The  motion  of  the  central  piston  in  the 
cavity  begins  when  the  shock,  driven  by  the  core  explo¬ 
sive,  reaches  the  front  surface  of  the  piston.  In  our 
calculations,  this  initial  motion  is  prescribed  based  on 
an  analysis  of  a  high  speed  framing  camera  record  of  the 
piston  motion  from  a  separate  experiment  (Reference  7). 
The  proper  momentum  of  the  piston  is  estimated  from  a 
l^-D  hydro  code  calculation  of  the  central  core  explosive 
event.  When  the  gas  arrives  in  the  cavity,  typically  1 
to  2  us  after  first  motion  of  the  piston,  the  motion  of 
the  piston  is  therafter  determined  by  the  difference  in 
pressure  across  the  piston.  The  initial  cavity  volume 
is  chosen  from  an  analysis  of  cavity  X-ray  data.  The 
cavity  walls  in  the  calculation  are  configured  to  give 
the  approximate  volume  as  a  function  of  time  determined 
from  cavity  X-ray  data.  Thus  the  energization  of  working 
gas  in  the  cavity  is  closely  linked  to  the  available  ex¬ 
perimental  data. 

With  the  initial  conditions  and  boundary  condi¬ 
tions  so  constrained,  the  mass  of  working  gas  is  the  re¬ 
maining  parameter  to  be  selected.  In  all  our  experiments 
we  begin  with  about  60  grams  of  gas  in  the  system.  From 
many  prior  tests  we  estimate  about  5%  of  this  is  lost  in 
the  driver  start-up  process.  From  cavity  X-rays,  we  es¬ 
timate  a  further  loss  in  the  region  where  the  gas  enters 
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the  end  cavity  (see  Figure  8  for  example) .  The  remaining 
gas--usually  25  to  35  grams — is  considered  available  for 
processing  in  the  cavity.  In  the  case  of  130-3  we  be¬ 
lieve  the  inlet  to  the  plasma  channel  collapses  rapidly 
cutting  off  the  flow  of  gas  to  the  channel  and  further 
reducing  the  mass  of  working  gas.  This  collapse  process 
is  modeled  by  using  a  collapse  trigger  level  and  radial 
collapse  history  based  on  the  2-D  code  calculations  of 
the  Voitenko  compressor  (Figure  28)  and  scaled  appro¬ 
priately  for  our  geometry. 
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Figure  28.  Calculation  of  the  Voitenko 
Generator  (Reference  8) 
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4.2.2  Comparison  of  Calculated  and  Measured 
Results 

The  calculated  and  measured  flow  velocity  and 
conductivity  histories  for  130-3  are  shown  in  Figures  29 
and  30.  The  calculated  velocity  is  higher  than  the  mea¬ 
sured  velocity  and  several  reasons  may  be  cited. 

1.  The  initial  entropy  state  of  the  source 
plasma  assumed  in  the  calculations  may 
be  too  high  since  the  plasma  is  pro¬ 
bably  processed  in  the  cavity  by  mul¬ 
tiple  shocks  rather  than  one  strong 
shock . 

2.  The  early  values  of  the  measured  ve¬ 
locity  may  be  lower  than  the  calculated 
1-D  velocity  because  of  the  non-unifor¬ 
mity  of  the  flow  and  the  manner  in  which 
the  velocity  gage  averages  the  flow. 

3.  At  later  times,  wall  expansion  could 
change  the  velocity  gage  electrode 
spacing  and  magnetic  field  thus  in¬ 
validating  the  gage  calibration. 
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Figure  29.  Calculated  and  Measured 
Flow  Velocity  Histories 
For  Shot  130-3 
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Fiaure  30  . 


Calculated  and  Measured  Conduct  i'/i 
Histories  For  Shots  130-3 
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4.  Boundary  layer  growth  and  flow  contam¬ 
ination  are  not  included  in  the  calcu¬ 
lations  and  if  important  would  tend  to 
reduce  the  calculated  velocity  at  late 
times . 

It  was  noted  during  the  course  of  several  cal¬ 
culations  that  the  calculated  velocity  was  insensitive 
within  limits  to  the  mass  of  plasma  initially  assumed  to 
be  available  for  processing  in  the  plasma  source  cavity. 
However,  lowering  the  initial  entropy  state  assumed  for 
the  calculations  would  directly  lower  the  calculated  ve¬ 
locity  since  the  flow  velocity  is  determined  primarily 
by  the  energy  density  of  the  reservoir  cas. 

The  calculated  form  of  the  conductivity  his¬ 
tories  and  decay  of  peak  conductivity  with  distance  rea¬ 
sonably  reflect  the  observed  behavior.  The  calculated 
plasma  pressure  and  internal  energy  histories  are  given 
in  Figures  31  and  32  and  can  be  used  to  find  the  other 
state  variables  in  the  1980  Livermore  equation-of-state 
for  argon.  In  the  calculation  14  grams  of  plasma  are 
driven  into  the  channel  before  the  inlet  is  completely 
collapsed.  The  effect  of  wall  expansion  can  be  aaged 
by  noting  that  the  total  enthalpy  of  the  plasma  pulse 
measured  at  any  point  along  the  diagnostic  channel  re- 

Wali  motion  is  too  slow  to 
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mair.s  constant  at  2.1  MJ. 


significantly  degrade  the  energy  of  the  pulse  since  it  is 
concentrated  at  the  high  velocity  leading  edge  of  the  pulse. 
The  compression  of  the  plasma  in  the  cavity  by  the  central 
piston  is  nearly  isentropic  to  a  peak  of  85  GPa.  Although 
this  compression  is  extremely  impulsive  (Figure  33) ,  the 
plasma  enters  the  cavity  at  a  high  entropy  level  from  the 
strong  shock  driver  process  and  shock  Mach  numbers  in  the 
cavity  are  never  high. 

The  calculated  and  measured  flow  velocity  and 
conductivity  histories  for  the  modified  plasma  source  of 
130-7  are  shown  in  Figures  34  and  35.  The  agreement  is 
fair  although  the  calculated  velocity  is  again  higher  for 
the  same  reasons  cited  previously.  In  the  calculation 
the  channel  inlet  is  not  forced  to  collapse  and  20  grams 
of  plasma  are  driven  into  the  channel  resulting  in  con¬ 
siderably  higher  plasma  densities  than  calculated  for  the 
130-3  experiment.  The  computed  flow  velocity  is  slightly 
lower  than  that  of  130-3  and  the  calculated  conductivity 
histories  are  in  general  slightly  higher  than  those  of 
130-3.  Peak  computed  cavity  pressure  is  25  GPa,  less  than 
one  third  of  that  calculated  for  130-3. 

The  X-ray  of  the  cavity  collapse  of  130-6  in¬ 
dicates  that  a  substantial  mass  of  plasma  is  trapped  by  a 
premature  collision  of  the  central  piston  with  the  cavity 
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Figure  34.  Calculated  and  Measured 
Flow  Velocity  Histories 
For  Shot  130-7 
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Figure  35.  Calculated  and  Measured  Conductivity 
Histories  For  Shot  130-7  _ 


walls.  In  the  calculation  it  was  necessary  to  reduce  the 
mass  of  plasma  entering  the  channel  to  about  8  grams  to 
achieve  reasonable  agreement  with  the  observed  flow  veloc¬ 
ity  and  conductivity  histories.  The  low  level  of  conduc¬ 
tivity  is  directly  related  to  the  low  average  gas  density 
in  the  channel. 

When  the  mass  of  plasma  in  the  channel  is  in 
rough  agreement  with  the  mass  expected  from  analysis  of 
the  cavity  X-rays,  the  agreement  between  calculation  and 
experiment  is  reasonable.  The  calculations  depend  strongly 
on  the  source  modeling,  on  the  validity  of  the  equation-of- 
state  and  on  the  applicability  of  the  Rogov  formulation  of 
conductivity.  Because  of  the  good  agreement  with  measured 
results,  we  are  confident  of  all  three  choices.  Channel 
t  ressure  history  measurements  are  planned  for  future  ex¬ 
periments  to  provide  additional  confirmation  of  our  hydro¬ 
dynamic  calculations. 

In  future  calculations  we  plan  to  incorporate 
ablative  boundary  layer  effects  into  the  description  of 
channel  flow.  As  this  effort  progresses  we  expect  to 
gain  additional  insights  into  the  hypothesis  that  early 
closure  of  a  blown  boundary  layer  can  have  a  limiting 
effect  on  the  duration  of  the  plasma  pulse.  We  are  also 
planning  to  add  an  MHD  power  extraction  model  to  the  code 
to  begin  to  determine  the  overall  effects  of  high  power 
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extraction  on  the  braking  of  the  plasma  flow.  We  also 
plan  to  review  the  choice  of  initial  plasma  state  assumed 
in  the  source  cavity.  We  are  considering  a  2-D  axisym- 
metric  hydro  calculation  to  more  accurately  assess  the 
processes  by  which  the  plasma  enters  the  cavity,  conver¬ 
ges  on-axis  and  flows  into  the  channel. 
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5.0  Conclusions  and  Recommendations 


5.1  Summary  of  Program  Accomplishments 

•  We  have  generated  and  measured  plasma  conditions 
of  extremely  high  energy,  electrical  conductivity 
and  flow  velocities.  Using  a  high  efficiency 
plasma  source  we  have  achieved  flow  velocities  of 
30  km/s  with  electrical  conductivities  as  high  as 
50  kS/m. 

•  An  excellent  analytical  framework  for  describing 
the  plasma  states  has  been  established.  The  1980 
Livermore  equation-of-state  for  ergon  is  now  used 
to  describe  the  thermodynamic  states  encountered 
in  the  plasma  source  and  channel  flow.  The  Rogov 
formulation  for  electrical  conductivity  is  used 
for  the  plasma  states  characteristic  of  the  chan¬ 
nel  flow. 

•  A  l^-D  hydrodynamic  code  technique  has  been  used 
to  model  the  plasma  energization  and  subsequent 
channel  flow.  Using  the  calculated  hydrodynamic 
flow  histories  and  the  Rogov  formulation  for  con¬ 
ductivity,  we  have  achieved  satisfactory  agree¬ 
ment  between  calculated  and  measured  conductivity 
histories  for  argon  plasmas. 
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•  In  the  channel  flow  plasma  states  characteris¬ 
tic  of  our  experiments,  the  range  of  measured 
electrical  conductivities  varies  from  50  kS/m 
at  peak  conditions  to  15  kS/m  at  the  end  of 
the  highly  expanded  plasma  pulse.  Xenon  exhi¬ 
bits  a  somewhat  higher  conductivity  than  argon 
or  air  at  the  same  thermodynamic  state. 

•  As  noted  in  previous  experiments,  measured  in¬ 
terelectrode  plasma  resistance  is  higher  than 
that  calculated  from  bulk  (search  coil)  conduc¬ 
tivity  measurements.  This  difference  is  tenta¬ 
tively  ascribed  to  differences  in  the  averaging 
process  by  the  two  measurement  techniques  in 

the  presence  of  a  less  conductive  boundary  layer. 

•  The  duration  of  the  plasma  pulse  is  relatively 
insensitive  to  source  size  and  source  geometry. 

We  hypothesize  that  pulse  duration  may  be  limited 
by  early  closure  of  a  blown  boundary  layer  in  the 
high  pressure  region  of  the  channel  near  the  plas¬ 
ma  source. 

•  Based  on  hydro  code  calculations  that  agree  with 
measured  flow  velocities  and  conductivities,  we 
estimate  a  peak  flow  enthalpy  of  a  terav/att  at 
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the  source  exit  and  0.2  terawatts  in  the  plasma 
diagnostic  section  of  the  channel.  The  total 
enthalpy  of  the  plasma  pulse  in  the  channel  is 
approximately  2.5  M J . 

•  We  have  measured  peak  electrical  power  densities 
of  3.6  TW/m'-T2  using  small  electrodes  in  a  field 
provided  by  permanent  magnets.  The  average  power 
density  over  the  first  25  us  of  the  plasma  is  2.4 
TW/nv--T~.  With  the  flow  conditions  already  achieved 
we  are  confident  of  extracting  electrical  power 
at  levels  of  tens  of  gigawatts  in  a  25  mm  channel 
with  an  externally  produced  10  tesla  field. 

A  summary  of  the  plasma  parameters  achieved  in  this 
ogram  is  given  in  Table  3. 
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Peak  Plasma  Source  Conditions 


Plasma  Pressure 


85  GPa 


Specific  Energy 


100  MJ/kg 


Flow  Power 

at  Channel  Entrance  -  1  TW 


Peak  Channel  Conditions  (at  20  diameters) 


Flow  Velocity 


30  km/s 


Conductivity 


50  kS/m 


Flow  Power 


0 . 2  TW 


Total  Pulse  Enthalpy  -  2.5  MJ 

(40  us  pulse) 

Peak  Electrical 

Power  Density  -  3.6  TW/m'-T- 

Average  Electrical 

Power  Density  -  2.4  TW/m'-T" 

(over  25  psec) 


Table  3.  Summary  of  Achieved  Performance  in 
Plasma  Source  ar.d  25  mm  Channel 
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5.2  Recommendations 


•  We  recommend  plasma  experiments  using  ar.  exter¬ 
nally  applied  magnetic  field  on  the  order  of  10 
tesla  to  achieve  high  electrical  power  output: 
and  to  investigate  the  moderate  plasma  flow  in¬ 
teractions  resulting  when  over  lOi  of  the  avail¬ 
able  hydrodynamic  power  is  extracted. 

•  Measurement  of  channel  flow  pressure  histories 
are  recommended  to  provide  further  confirmation 
of  channel  plasma  conditions. 

•  Investigations  should  be  continued  on  the  limit¬ 
ing  effect  of  a  blown  boundary  layer  on  the  pulse 
duration  of  plasma  flow  in  the  channel.  This 
analysis  may  also  help  clarify  the  effect  of 
boundary  layer  growth  on  interelectrode  resis¬ 
tance  . 

•  Further  development  of  the  plasma  source  should 
be  undertaken  to  increase  the  mass  of  plasma 
driven  into  the  channel  since  increased  flow- 
densities  directly  increase  the  level  of  plasma 
conductivity . 
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•  Continued  development  of  the  1H-D  computations 
are  recommended  to  model  the  overall  effects  of 
MHD  power  extraction  on  the  plasma  flow  and  its 
thernodynami c  s ta te . 

•  A  program  to  study  the  detailed  3-D  interaction 
of  the  plasma  flow,  magnetic  field  and  electric 
power  extraction  circuit  should  be  initiated  to 
form  the  basis  of  a  fully  interactive  3-D  MHD 
computer  code. 
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